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Abstract 

In this paper the description of the dielectric and ferroelectric properties of BaTiOs single crystals 
using Landau phenomenological thermodynamic potential was addressed. Our results suggest that 
when using the sixth-power free energy expansion of the thermodynamic potential, remarkably 
different values of the fourth-power coefficient, (3 (the coefficient of terms), are required to 
adequately reproduce the nonlinear dielectric behavior of the paraelectric phase and the electric 
field induced ferroelectric phase (£^ft||[001]), respectively. In contrast, the eighth-power expansion 
with a common set of coefficients enables a good description for both phases at the same time. 
These features, together with the data available in literature, strongly attest to the necessity of 
the eighth-power terms in Landau thermodynamic potential of BaTiOs. In addition, the fourth- 
power coefficients, (3 and ^ (the coefficient of PfPj terms), were evaluated from the nonlinear 
dielectric responses along [001], [Oil], and [111] orientations in the paraelectric phase. Appreciable 
temperature dependence was evidenced for both coefficients above Tc- Further analysis on the 
linear dielectric response of the single domain crystal in the tetragonal phase demonstrated that 
temperature dependent anharmonic coefficients are also necessary for an adequate description of 
the dielectric behavior in the ferroelectric phase below Tc- As a consequence, an eighth-power 
thermodynamic potential, with some of the anharmonic coefficients being temperature dependent, 
was proposed and compared with the existing potentials. In general the potential proposed in 
this work exhibits a higher quality in reproducing the dielectric and ferroelectric properties of this 
prototypic ferroelectric substance. 
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I. INTRODUCTION 



Structural phase transformation phenomena are an important issue in sohd state physics. 
Both phenomenological and microscopical approaches have been playing comparable roles 
in understanding the phase transition nature of the magnetics, superconductors, and ferro- 
electrics. Particularly for ferroelectrics, considerable attention is paid to the phenomeno- 
logical Landau-Devonshire theory in parallel with advancing techniques of first-principle 
calculations. In the phenomenology of ferroelectrics, the issue recently attracting atten- 
tion was the role of high-order anharmonic polarization terms in the Landau free-energy 
expansion.-!^ In these works, the terms containing eighth (or higher) power of polarization, 
traditionally neglected in the phenomenological framework, was demonstrated to be vitally 
important for the description of the phase diagrams of perovskite ferroelectrics (e.g. for the 
recently discovered monoclinic phase of some complex perovskites in the vicinity of the mor- 
photropic phase boundary).'^ The physics behind this effect is the need of high-order terms 
for adequate description of the symmetry of the problem. Another demarche in the field has 
been recently undertaken by Li et ai, who have revisited the phenomenological theory of 
classical ferroelectric BaTiOs.- These authors suggested incorporation of the eighth power 
terms in the Landau expansion of BaTiOs in order to eliminate the temperature dependence 
of the anharmonic coefficients and to adapt the description to the situation of large compres- 
sive strains. In contrast to Refs. |1|J3, here no symmetry arguments have been involved. The 
message by Li et al. was that the Landau expansion containing eighth power terms (with 
temperature independent anharmonic coefficients) is as efficient in the description of many 
properties of BaTiOa as the traditional sixth order expansion with the temperature depen- 
dent anharmonic terms.-i^ This conjecture is of the interest since the strong temperature 
dependence of the anharmonic coefficients in the BaTiOs Landau expansion is in conflict 
with the displacive nature of ferroelectricity in the material, at least under the common 
assumption that the critical fluctuations are weak in this system. ^'^ Interestingly, in their 
analysis of the Landau potential using first-principles calculations, Iniguez et al. obtained 
that the sixth order expansion accounts for the main features of BaTiOs phase diagram, but 
that all coefficients in the expansion have nontrivial temperature dependences.- 



3 



All in all, these results pose two questions. First, are higher than sixth power terms in the 
BaTiOs Landau expansion really needed for the correct physical description and are not just 
a matter of convenience? Second, does the incorporation of such terms enable elimination 
of the strong temperature dependence of the anharmonic coefficients in the BaTiOa Lan- 
dau expansion? Addressing these questions, the present paper demonstrates by analyzing 
the dielectric response that, both the eighth-power terms and the temperature dependence 
of the anharmonic coefficients (at least of the quartic coefficients) are essential for the ad- 
equate description of the thermodynamic behavior of BaTiOa, especially on crossing the 
cubic-tetragonal phase transition. Including these necessities, an improved Landau phe- 
nomenological potential is proposed and compared with the existing models. Preliminary 
results on this subject have been published elsewhere.— 

II. IMPORTANCE OF THE EIGHTH-POWER TERMS 

A. The effect of the eighth-power terms on dielectric nonlinearity 

Besides the phase stability, the nonlinear dielectric behavior is another characteristic 
that directly evidences the contributions of the anharmonic terms in Landau thermodynamic 
potential. To demonstrate the crucial role of the eighth-power term in the Landau expansion 
we firstly address a simple case, namely the dielectric nonlinearity of BaTiOa with respect 
to the electric field E applied along the [001] axis, in the vicinity of the cubic/tetragonal 
phase transition temperature, Tc- In this case only the [001] component of the polarization, 
P = P3, is involved, so that the Landau expansion of the Gibbs potential containing the 
eighth-power term reads 




(1) 



Accordingly we have for the equation of state and c-axis permittivity ec- 



E = aP + pP^ + -fP^ + 6P 



(2) 



-1 



a + 3pP^ + 5-fP^ + 75P 



.6 



(3) 
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These equations provide a description of the non-hnear dielectric response of BaTiOs in 
terms of the eighth-power expansion. One can compare this description with that in terms 
of the sixth-power expansion 

E = aP + + (4) 

e-^ =a + 3/3' + 5^'P^ (5) 

which have been traditionally used in the field. These descriptions are clearly different. 
Formally, as one can readily check, the situation described by the expansion Eq. ([2]) and 
([3]) with the polarization-independent coefficients /3, 7, and 6 corresponds to polarization- 
dependent coefficients /5' and 7' in Eq. (jlj) and ([5]): 

f3' = (3- 5P\ 7' = 7 + 25P2 (6) 

In general, the consideration of the polarization-dependent anharmonic coefficients in the 
Landau expansion has no physical sense. However, if a not too wide interval of polarization 
variation is experimentally addressed with a finite measurement accuracy, coefficients 13' 
and 7' defined by Eq. ([6]) can be considered as polarization-independent when fitting the 
dielectric non-linearity. At the same time, serious problems may arise in situations when a 
large variation of polarization is involved, e.g. the "jump" of the polarization on crossing 
the first order phase transition. To verify this point we consider the case just above Tc, 
where the ferroelectric phase can be induced by an electrical field. As is schematized by 
the ^ segment of the "double hysteresis loop" in Fig. [H the polarization "jumps" 
at the critical field Ec-, whereas below or above this value the polarization varies slowly 
with bias field. At electric fields corresponding to the paraelectric phase, the polarization 
varies from zero to the critical value Pc- Consequently the contributions of the polarization 
dependent terms can be estimated from Eq. ([6]) to be not more than 5Pq and 25P^ to (3' 
and 7', respectively. In contrast, at electric fields corresponding to the ferroelectric phase, 
the polarization ranges in a much higher interval of not less than Pq? the minimum value of 
polarization in ferroelectric phase. Accordingly the polarization dependent correction terms 
in Eq. ([6]) will be expected to be not less than 5Pq and 25Pq , respectively. In Table [U 
we list the numerical values of the polarization dependent terms in both phases and the 
typical values of the corresponding coefficients. The value of 5 follows Li et al.'s suggestion.- 
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Apparently the values of the polarization dependent terms strongly depend on the interval 
over which the polarization varies. In ferroelectric phase, the values of the "correction" 
terms are comparable to the typical values of the corresponding coefficients, whereas they 
are negligible in paraelectric phase. 

Thus from the above analysis, we conclude that if the coefficients of the eighth-power 
terms have a value of the order of that suggested by Li et al, the anharmonic polarization 
coefficients /?' and 7' of BaTiOs estimated in terms of the sixth-power expansion should 
exhibit appreciable jumps on crossing the field induced phase transition. Such "phase sensi- 
tivity" of the nonlinear coefficients of the sixth-power expansion enables the detection of the 
eighth-power terms, suggesting a simple and direct way to assess the role of the eighth-power 
terms. 

Keeping this in mind, we studied the nonlinear dielectric behavior of BaTiOs single 
crystals in paraelectric phase and in electric field induced ferroelectric phase (-E';,//[001]) 
at a temperature a few degrees higher than Tc- The experimental details can be found 
in our previous work.— Figure [2] shows the experimental data measured at 135°C (for this 
crystal, the phase transition temperature Tc = 130°C). Little hysteresis effect related to the 
relaxation phenomena is found within either phase, indicating the high insulating quality 
of the crystal and the space charge free state of the crystal-electrode interface.— The sixth- 
order expansions, Eqs. (jlj) and are used to fit the curves with the common value of a, 
and remarkably different (3 are obtained for either phase, as listed in Table [ITl As for 7, 
we cannot determine it exclusively from the nonlinear dielectric behavior of the paraelectric 
phase. The position of Ec, where the ferroelectric phase is induced by the bias field, is 
used to obtain an estimate, which differs considerably from that derived in the ferroelectric 
phase. The ferroelectric-phase coefficients have a serious problem in describing the dielectric 
properties of the paraelectric phase, and vice versa. As illustrated in the inset of Fig. [21 the 
eighth order expansion with a common set of coefficients enables a good fit in both phases 
at the same time. As-estimated values of /3, 7, and 6, agree well with the set suggested by 
Li et alA 

Furthermore this approach is extended to a wider temperature interval around Tc- As 
shown in Fig. [31 the quartic coefficient determined with the sixth-power expansion exhibits 
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strong phase dependence, characterized by a "jump" on crossing the phase transition no 
matter whether this transition is induced by electric field or temperature. In contrast the 
values determined with the eighth-power expansion are phase independent as Landau theory 
coefficients must behave. The "phase sensitivity" of (3 determined by the sixth-power ex- 
pansion, and its elimination by the incorporation of the eighth-power term, strongly attest 
to the essentiality of the eighth-power terms in BaTiOs Landau expansion. 

B. Literature available evidence 

In the preceding section, we show that the eighth-power terms are essential for the ad- 
equate description of the dielectric nonlinearity of BaTiOa, especially in the ferroelectric 
region in the vicinity of the cubic-tetragonal phase transition temperature. Interestingly, 
this point was not realized for a long time. The probable reason was that the anharmonic 
coefficients were in most cases determined by dealing with the properties (e.g., dielectric 
permittivity,—'^ birefringence,— and phase instability^^) of the paraelectric phase using the 
sixth-power expansion, whereas the dielectric behavior or the polarization of the ferroelec- 
tric phase were not considered. In these cases the contribution of the eighth-power terms is 
negligible and cannot be detected. Values of the coefficients reported in these early studies 
are actually equivalent to our nonlinear coefficients determined from the paraelectric phase 
only, as explained above. When the properties of the ferroelectric phase, particularly the 
very high polarization are involved in determining the anharmonic coefficients of the sixth- 
power expansion, quite different values for the fourth and sixth order coefficients can be 
obtained according to our preceding discussion. Actually this divergency can be found in 
Gonzalo et a/.'s work,^^ where the critical polarization values, at which the paraelectric or 
the ferroelectric phase loose their stability, have been employed in evaluating the anhar- 
monic coefficients. The values of their coefficients are nearly the same as our coefficients 
determined from the ferroelectric phase. In Table IIIII we collected values of the coefficients 
available from the literature, which were determined in the vicinity of Tc- The method 
of determining coefficients is brieffy described for each case. Appreciable difference can be 
found when comparing Gonzalo et a/.'s data and those derived from the properties of the 
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paraelectric phase. Unfortunately this difference of ~ 100% in magnitude for (3 and ~ 300% 
for 7 were ignored by the authors though it provides a strong signal that the eighth-power 
terms cannot be neglected. 

Additional evidences for the necessity of the eighth-power terms can be found from the 
linear dielectric responses on crossing the cubic-tetragonal phase transition. At the first 
order phase transition temperature Tc, under zero bias field, we have 

AG(P5) = \aPl + \pPl + \lPl + \^Pl = (7) 

E{Ps) = aPs + PPI + P| + 5P] = (8) 

e-i = a + 3(3Pl + 5Pl + 76Pl (9) 

= « = «o(rc - To) (10) 

where Ps is the spontaneous polarization of the tetragonal phase at Tq. ec,c and ec,T 
represent the c-axis dielectric permittivity of the cubic and the tetragonal phase on crossing 
the phase transition, respectively. ao=l/(Ceo) is the temperature derivative of a and C, Tq 
the Curie- Weiss constant and Curie- Weiss temperature, respectively. 

Accordingly, two quantities, namely the variation of the dielectric permittivity and the 
variation of the temperature derivative of the inverse permittivity along c-axis on crossing 
the ferroelectric phase transition, can be analytically derived as 

i? = !£i£=4 + ^ (11) 
ec,T a 

^'--^^-8 + ^^-3^^-5^^-... (12) 

c^(l/ec,c) 47 + 95P| aodT dT ^ ' 

dT 

Using the coefficients values and their temperature derivatives of the existing Landau 
potentials,-'^'^ we can estimate the numerical values for the two quantities. Table [IV] lists 
the numerical results calculated by the analytical derivations with sixth- and eighth-power 
expansions. In case of the sixth-power expansion, namely 5 = 0, P = 4 and R' ~ 8. In 
contrast, the existence of 5 that has a value of the order of that suggested by Li et ai, 
leads to appreciable larger values for both quantities. Available experimental data are also 
collected for comparison. Apparently most of the experimental values are considerable larger 
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than those derived by the sixth-power expansion. Considering the possible experimental 
inaccuracy in determining the lattice permittivity in the ferroelectric phase due to the not 
fully eliminated domain contribution, the real values of R is expected to be even larger.— 
The deviations of experimentally determined R and R' from those calculated from the sixth- 
power expansion, affords additional evidence for the essentiality of the eighth-power term. 

In summary the necessity of the eighth-power terms in BaTiOs Landau expansion is 
demonstrated by our dielectric nonlinearity measurements and literature available data. 
Keeping in mind that the existence of the eighth-power terms leads to the polarization 
dependence of the anharmonic coefficients in the sixth-power expansion according to Eq. ([6]), 
the temperature dependence of the quartic coefficients, which was usually determined with 
the sixth-power expansion,— might be logically ascribed to the eighth-power terms as the 
polarization around Tc is strongly temperature dependent. In next section we will address 
the temperature dependence of the anharmonic coefficients in the presence of the eighth- 
power terms and show that the incorporation of the eighth-power terms cannot eliminate 
the temperature dependence of the fourth order anharmonic coefficients. 

III. TEMPERATURE DEPENDENCE OF THE ANHARMONIC COEFFICIENTS 
A. Temperature dependence of the quartic coefficients 

Though the existence of the non-zero eighth-power terms leads to polarization dependent 
corrections to the quartic and sextic coefficients, however, these corrections are negligible in 
the paraelectric phase regime, where the contribution of the anharmonic terms is dominated 
by the quartic terms. Consequently the values determined with the sixth-power expansion 
of paraelectric properties have a high accuracy. As can be seen from Fig. [3], the quartic 
coefficient determined by the sixth-power and eighth-power expansion agree well when only 
paraelectric data are considered. Therefore we cannot disregard the temperature depen- 
dence of (3 suggested by Drougard^^ and Meyerhofer— according to the nonlinear dielectric 
or birefringence behavior of the paraelectric phase. A linear temperature dependence was in- 
dependently revealed for /3 at temperatures above but not far above Tc- Figure ?? illustrates 
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the quartic coefficient /3 as a function of temperature in the vicinity of Tc- An apparent 
temperature dependence can be found in the paraelectric phase, well agreeing with Drougard 
and Meyerhofer's results. Linearly extrapolating to higher temperatures shows that f3 van- 
ishes at around 185°C, which, however, is inconsistent with our previous observations of 
around 165°C.^° The difference in the zero cross-over temperature of f3 may be influenced 
by many factors, e.g. the conductivity, charge accumulating at the oxide-electrode interface, 
and processing history of the sample, all of which can strongly affect the accuracy of the 
nonlinear measurements. It is important, however, that both our previous and present work 
show qualitatively the same behavior of /3, which is temperature dependent and changes sign 
at some temperature well above Tc- In the ferroelectric phase, the f3 values determined with 
the eighth-power expansion and temperature independent 7 and 6 exhibit a continuous and 
similar temperature dependence as in paraelectric phase (see Fig. [3]). However this tempera- 
ture dependence is, by itself, not convincing since we cannot extract sufficient information to 
verify the temperature dependence of the higher-than-fourth power anharmonic coefficients 
merely from the dielectric nonlinearity measurements. The necessity of the temperature 
dependent anharmonic coefficients will be demonstrated in the next subsection. 

The value of another quartic coefficient, ^ [see Eq. (fT3l) ]. can be evaluated by combining 
the nonlinear dielectric responses along [001], [Oil], and [111] orientations of the paraelectric 
phase. In this case the anharmonic contributions are dominated by the quartic terms, and 
the higher-than-fourth-power terms can be omitted.— The free energy expression can be 
simplified as 

AG = la{P^ + P| + P|) + \f3{P,' + + Pi) 

+\^{P!Pi + PiPi + PiPl) - EiPi - E2P2 - E3P3 (13) 

where Pi, P2, P3 and Ei, E2, E^ are cartesian components of polarization and electric field, 
respectively. The nonlinear dielectric response along three orientations can be deduced by 
successive differentiation:— 

Ehki = aPhki + f^hkiPhki (14) 
e,,\ = a + 3PHkiPLi (15) 
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where the effective nonhnear coefficients (3hki are hnear combinations of [3 and ^: 



/^ooi — 

Poll = liP + (16) 
All = + 

Figure ?? illustrates the values of ^ as a function of temperature. The values determined 
from various combinations agree with each other with an accuracy of 10% and show consis- 
tent temperature dependence over the studied temperature interval. However, our results 
determined from dielectric nonlinearity are quite different from those proposed in Ref. Il4l . 
where the values of ^ were directly estimated by fitting the electric field dependence of the 
induced polarization of the paraelectric phase. Both the value itself and its temperature 
derivative have opposite signs in comparison with those reported by Meyerhofer— in the 
vicinity of Tq. We believe that our positive ^ values are correct for the following reasons. 
Both [Oil] and [111] oriented crystals exhibit a decreasing electric field dependence of the 
dielectric permittivity, suggesting positive values of the effective nonlinear coefficients, /?oii 
and Pui. In contrast, a negative ^ definitely results in negative /3oii and/^m, which is 



opposite to our experimental data. A negative value of ^ suggested in Ref. [ij is actually in 
contradiction with the birefringence measurements reported in the same paper.— Finally, ^ 
must exhibit a value larger than the absolute value of f3 to fit the positive /?oii- A smaller 
value of C,, as was implicitly assessed by fitting the tetragonal/orthorhombic phase transition 
temperature with the sixth-power expansion,-i^i2^ cannot describe the dielectric nonlinearity 
of [Oil] crystals in the paraelectric phase. 

Therefore, we show that even in the presence of eighth-power terms, which are actually 
not important for the paraelectric phase, both of the quartic coefficients definitely have to be 
temperature dependent, at least above Tq. The temperature dependence of the higher-power 
coefficients, however, cannot be precisely determined merely from dielectric nonlinearity 
results. More measurements, e.g. the polarization as a function of electric field, are needed 
to clarify this point. 
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B. Temperature dependence of the anharmonic coefficients in tetragonal phase 



In ferroelectric phases below Tc, the contributions of the terms higher than quartic are 
comparable to those of the quartic terms.— Neither the coefficients themselves nor their tem- 
perature dependence can be determined from dielectric nonlinearity in a reasonable accuracy. 
Here we propose a simple way to cross-check the necessity of the temperature dependent 
anharmonic coefficients in tetragonal phase. Firstly we assume that the Landau coefficients 
are temperature independent except the quadratic coefficient, a, which is linearly tempera- 
ture dependent following Curie- Weiss law. The equation of state and c-axis permittivity ec 
of the tetragonal phase can be simplified as: 

E = ao{T-To)P + g{P) (17) 

e-^ = ao{T-To)+g\P) (18) 

where g{P) summarizes all the contributions from the higher order anharmonic terms. Under 
zero electric field, = 0, P = Pg- Assuming that the only temperature dependence in g 
comes from P5, one can readily derive 

(\W d P 

— = = aoP5 + [«o(T - To) + g\Ps)]^ (19) 



Thus we have 



Combining Eq. (fTS!) and (!20|) . we get 



1 dlnPs ^dlnPs , . 



aoCo dT dT 

where e^^c represents the relative dielectric permittivity along the direction of the spontaneous 
polarization P5 in tetragonal phase. Thus by comparing the quantity — e^^c/ {dlnPs/dT) with 
C, we can readily verify the validity of the temperature independent assumption. 

In Fig. [5] we compare the two quantities as a function of temperature in the tetragonal 
phase. The experimental data are extracted from Ref. 
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where the measurements were 



performed on identical single domain BaTiOs crystals. Appreciable difference (~ 50% of 
C) can be found between these two quantities, attesting to the invalidity of the tempera- 
ture independent assumption. Namely, in ferroelectric phase, the temperature dependent 
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anharmonic coefficients are indeed necessary for an adequate description of the thermody- 
namics of BaTiOs- However, detailed information about the temperature dependence of 
each coefficient cannot be obtained by this simple method. 



IV. IMPROVED THERMODYNAMIC POTENTIAL FOR BATIO3 SINGLE 
CRYSTALS 

In the preceding sections we demonstrate that not only the eighth-power terms, but 
also the temperature dependent anharmonic coefficients (at least temperature dependent 
quartic coefficients), are essential in BaTiOs Landau thermodynamic potential. Hereinafter 
we propose the following improved potential for this prototype substance by introducing 
temperature dependent coefficients into the eighth-order expansion. For convenience we use 
the same notations as in Ref. 

AG = ai(p2 + p| + Pl) + an (Pi' + + P^) 

+au{P^P^ + P^Pi + P^P^) + ainiP^ + P^ + P^ 

+au2 [Pm + P3) + PiiPl + Pt) + ^3 (^1' + P2 



+ai23Pi'P|^3 + + P! + P'. 



3 , 



+«iii2 [p^Pi + Pi) + p|(p| + Pi) + Pl{P! + Pl 

+a,,,,{PtP^ + P^Pt + PtPt) 

+«ii23(Pi'P|^3 + p!p2pI + pMpI) 

-E,P, - E2P2 - E3P3 (22) 

where the coefficients ai, an, ai2, am, and aim corresponding to the preceding notations 
^a, -2^, ^7, and -^5, respectively. The temperature dependent ai, an, and ai2, are 
directly determined from our linear and nonlinear dielectric measurements above Tc- The 
linear dependences are simply extrapolated down to low temperatures. The temperature 
dependence of am is estimated by ffiting the spontaneous polarization of the tetragonal 
phase, with a temperature independent aim that is evaluated from the dielectric nonlinear- 
ity above Tc- A linear dependence is revealed for am and extrapolated to low temperatures. 
The rest of the coefficients are assumed to be temperature independent. aii2 and aiii2 are 
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estimated by fitting the dielectric permittivity perpendicular to the spontaneous polariza- 
tion, ea, as a function of temperature in the tetragonal phase. 0:1122 is evaluated to fulfill 
the tetragonal/orthorhombic phase transition temperature. 0123 and 01123 are adjusted to 
fit the orthorhombic/rhombohedral phase transition temperature and its field dependence.— 
The numerical values of the nonlinear coefficients described above are listed in Tab. |Vl For 
comparison those of the two existing potentials proposed by Bell^ and Li et alA are also 
presented. 

In the first glance, our model is just a combination of those proposed by Li et al. and 
Bell.->^ However it should be kept in mind that (i) the eighth-power terms used in Li et 
a/.'s potential are not just a convenience but an essentiality in eliminating the "phase sen- 
sitivity" of the anharmonic coefficients of the sixth-power expansion; (ii) even with the 
eighth-power expansion, the temperature dependence of the anharmonic coefficients can- 
not be disregarded, as has been generally demonstrated by the dielectric properties of the 
tetragonal phase (Fig. [5]), and directly revealed by our dielectric measurements for the two 
quartic coefficients (Fig. H]). Disregarding these temperature dependence of the quartic 
coefficients and the eighth order terms leads to an inadequate description of the dielectric 
behavior, even in a qualitative way. In contrast, incorporating these changes, a high quality 
in reproducing the dielectric behavior of BaTi03 single crystals is achieved, particularly in 
the vicinity of the cubic-tetragonal phase transition. Firstly, above Tq our model is capable 
in quantitatively fitting the dielectric nonlinearity as a function of temperature and crys- 
tallographic orientation (see Fig. [2]and|l]), which in principle cannot be well described by 
either of the other potentials. Secondly, in the tetragonal phase the linear dielectric con- 
stant and the spontaneous polarization can be precisely reproduced. Figure El compares the 
experimentally measured and the theoretically calculated 1 / e^r- as a function of temperature 
around the paraelectric/ferroelectric phase transition. In spite of the great difference exist- 
ing among the investigated crystals (e.g. the Tc varies from 108 °C for Merz's crystal^ to 
130 °C of our crystal), the measured dielectric permittivity concide on departing from Tc- 
The curve calculated with the potential proposed in this work agrees well with the experi- 
mental data, whereas those calculated with Li et a/'s and Bell's coefficients show remarkable 
divergence with the experimental values in tetragonal phase. Note that when determining 
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the anharmonic coefficients we have never used the experimental data measured below 
Tc. 

The low temperature ferroelectric properties are also calculated with the three potentials 
listed in Tab |Vl The three consecutive phase transition temperatures are yielded with the 
current potential as 404 284 and 189 respectively. In Fig. [TJ we present the 
spontaneous polarization along pseudo cubic [001] of BaTiOs as a function of temperature. 
Some of the experimental data available in literature are also presented for comparison. In 
the tetragonal phase, the three potentials have similar quality in reproducing the sponta- 
neous polarization, both in its absolute value and in its temperature dependence. In the 
orthorhombic phase none of these models can well describe the temperature dependence of 
the spontaneous polarization. The existing potentials (Li et aVs and Bell's) give better es- 
timations of the average value of the spontaneous polarization. In the rhombohedral phase, 
however, our potential exhibits a better accuracy than the other two potentials. Again we 
emphasize that the experimentally measured polarization values of the rhombohedral phase 
have not been used in determining the anharmonic coefficients. As for the dielectric prop- 
erties of the low temperature ferroelectric phases, experimental data available in literature 
are rather scattered and sometimes are not self-consistent.— "^^i^^ Tentatively we measured 
a poled [111] crystal at cryogenic temperatures. At 20 K the relative dielectric constant is 
about 40. Considering the possible incompletely eliminated domain contributions, the real 
value seems to be even less. The corresponding calculations give 30, 81, and 68 when using 
our, Li et a/.'s and Bell's potentials, respectively. Apparently the thermodynamic potential 
proposed in this work exhibits a better accuracy. 



V. CONCLUSIONS 



In conclusion, we demonstrate by analyzing the dielectric nonlinearity that both the 
eighth-power terms and the temperature dependent anharmonic coefficients are necessary in 
BaTiOa Landau phenomenological potential for the adequate description of the thermody- 
namic behavior of this substance. Accordingly, an improved Landau potential, characterized 
by an eighth-order Landau expansion with some temperature dependent anharmonic coef- 
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ficients, is proposed and compared with the existing potentials. In general the proposed 
potential exhibits a higher quality in reproducing the dielectric and ferroelectric properties 
of this prototype ferroelectric substance. 
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TABLE I: Estimated values of the polarization dependent terms in Eq. ([6]) and typical values of 
the corresponding coefficients. The value of 6 is taken from^ 



Quantities Paraelectric Ferroelectric Units 

P ^ 0.04 ^ 0.16 (C/m^) 

5P^ ^ 5 X 10"^ ^ 2 X 10® {Vm^C~^) 

2dP^ ^ 2 X 10^ ^ 1.5 X 10^° (Vm^C-^) 

-P (6 ~ 9) X 10^ (Vm^C-^) 

7 (1 - 2) X 10^° {Vm^C-^) 
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TABLE II: Anharmonic coefficients determined from the field dependence of dielectric permittivity 
shown in Fig. [2j The values proposed by Li et al. are also listed for comparison. 



Order Phase(s)/Author -(3{\w.^C-^) -i{\m^C-^) 5{\m^^C''^) 

6 Paraelectric 8.6 x 10*^ 1 x 10^° 

6 Ferroelectric 18.7 x 10^ 8.2 x 10^° 

8 Both 8.6 X 10*^ 0.8 x 10^° 3.4 x 10^^ 

8 Li et alA 8.4 x 10*^ 0.77 x 10^° 3.0 x 10^^ 



"The value of 7 in paraelectric phase cannot be determined only by fitting the electric field dependence of 
the permittivity. The present value is estimated by combining the dielectric nonlinearity and the position 
of the critical field Ec- 
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TABLE III: Literature available values of the anharmonic coefficients for BaTiOa determined with 
the sixth-power expansion in the vicinity of Tc- 

Authors —f3(Vm^'C^^) 7{Vm'*C^^) Principles of the determination 

5.4 X 10* 1.7 X 10^'^ Instability of the paraelectric phase under bias field 

(4 ~ 8) X 10* NA Polarization dependence of the dielectric permittivity 

in paraelectric phase 

(7 ~ 10) X 10* (2 2.5) X 10^^ Field dependence of the polarization in paraelectric 

phase 

6 X 10* NA Field dependence of the dielectric permittivity in para- 

electric phase 

(15 ~ 20) X 10* (5 ~ 10) X 10^^ Instability of the paraelectric and ferroelectric phase 

above Tc" 

"In this work, the critical polarization (or the polarization at the destabilizing electric field) of either phases, 
were determined from double hysteresis loops. These values in practice strongly depend on the measuring 
and sample conditions, e.g. the switching frequency and electrode/crystal interface natures. Even though, 
the differences are appreciable when comparing with those determined from paraelectric properties. 



MerJi^ 
Drougardi^ 

Meyerhofei™ 

Kaczmareki^ 

Gonzala^ 
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TABLE IV: Analytical calculated and experimentally determined values of R and R' for BaTi03 
single crystals. The anharmonic coefficients and their temperature derivatives are taken from 
literature^i^!^. The spontaneous polarization are taken from literatur©i£. 



Quantities 


6th-power 


8th-power 


Experimental data 




expansion 


expansion 




R 


4 


5.3 


6.2^ 7.3^ 3.3^^, 6.1^ 5.2/ 


R' 


7.9^ 


9.3 


125, 9.4/ 



°As can be seen that the contributions by the temperature dependence of the anharmonic coefficients to 

R' are quite weak and in an opposite direction in comparison to that by the eighth-power term. 
^RefJ^ 
'^Ref^ 
'^Ref-H 

^This work 
SRef^ 
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TABLE V: Coefficients of BaTiOa Landau thermodynamic potential in Eq. (j22p . where T 
temperature in K. 



CoefBcients 


This work 




Li et al. " 


Bell and Cross 


h 


Units 




3.61 X 105(T 


- 391) 


4.124 X 105(T - 388) 


3.34 X 10^ (T - 


381) 




"11 


-1.83 X 10^ - 


h4 X lO^T 


-2.097 X 10» 


-2.045 X 10^ 4 


- 4.69 X 10<3T 


VmJ'C-^ 


Ol2 


-2.24 X 10^ - 


h 6.7 X lO'^T 


7.974 X 10® 


3.23 X 10** 






"111 


1.39 X 10^0 - 


3.2 X lO^T 


1.294 X 10^ 


2.445 X 10^" - 


5.52 X lO^T 


Vm^C'-^ 


"112 


-2.2 X 10** 




-1.950 X 10^ 


4.47 X 10^ 






"123 


5.51 X 10^° 




-2.500 X 10^ 


4.91 X 10^ 




Vm^C'-^ 


"1111 


4.84 X 10^0 




3.863 X 10^0 









"1112 


2.53 X 10^^ 




2.529 X 10^° 









01122 


2.80 X 10" 




1.637 X 10^° 









Q1123 


9.35 X 10^" 




1.367 X 10^" 
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FIG. 1: Schematic of the phase transition induced by an electrical field applied along [001] axis 
above Tq for BaTiOs single crystal. In paraelectric phase the polarization varies with field below 
Pc, while in ferroelectric phase the polarization sits in the interval above Pq- 



FIG. 2: (Color online) Field dependence of the dielectric permittivity for BaTiOa single crystal at 
135 °C. The open squares represent the measured values, and the dotted line (blue) is to guide 
the eyes. Remarkably different /? are obtained for either phase when fitting with Eq. (HD and (0), 
as listed in Table [H The dashed (red) and solid (green) lines are calculated with the coefficients 
derived from the paraelectric and ferroelectric data, respectively. The ferroelectric coefficients have 
a serious problem in describing the dielectric properties of the paraelectric phase, and vice versa. 
The inset shows the much better fit quality (black dashed-dotted line) of an eighth-order expansion 
for both phases. 



FIG. 3: (Color online) The quartic coefficient determined with sixth and eighth-power expansions 
in the vicinity of Tc- Phase dependent (3' values are derived when using sixth-power expansion, 
whereas phase independent /3 can be obtained with the eighth-power expansion. The /3 values 
below Tc are estimated by fitting the dielectric nonlinearity with constant 7 and S. However this 
is less convincing since the temperature dependence of 7 and 6 have been arbitrarily ommitted. 
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FIG. 4: (Color online) Temperature dependence of the quartic coefficients (3 (a) and ^ (b). The 
values of ^ are determined by combining the dielectric nonlinearity measurements along [001], 
[Oil] and [111] orientations, and the values and temperature derivatives calculated from three 
combinations agree well with each other. 



FIG. 5: (Color online) Divergency between the quantity — ^^^^p '^^j- ^'^d ^ ™- ferroelectric 
tetragonal phase. The values of dielectric permittivity and spontaneous polarization are extracted 
from literaturei^. 



FIG. 6: (Color online) Inverse dielectric permittivity along the spontaneous polarization in the 
vicinity of T^, a comparison between the calculations and the experimental data. For convenience 
the phase transition temperatures are normalized to 120°C. The scattered dark marks are ex- 
perimentally measured by different researchers on different crystals: "+" by Drougard^, "Q" 
by Merai^, and "□" in this work. The dashed (blue), dash-dotted (green), and solid (red) lines 
are calculated with Li et aZ.'s, Bell's and the newly-proposed potential in this work, respectively. 
Note that despite the great difference in the source sample, the permittivity exhibits tremendous 
consistence on departing from Tc- The new potential has a much better quality in fitting the 
experimental data than the two existing potentials do. 
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FIG. 7: (Color online) Spontaneous polarization along pseudo cubic [001] of BaTiOs, a comparison 
between the calculations and the experimental data. The scattered marks represent the experi- 
mental data: and "A" are extracted from Iiteraturei^i2^i2^, respectively. The dashed 
(dark), dash-dotted (red) and solid (blue) lines are calculated with Li et aZ.'s, Bell's and the newly- 
proposed potential in this work, respectively. The polarization values of the rhombohedral phase 
are exactly reproduced by the improved thermodynamic potential. 
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